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Transient Hamiltonian chaos in the cavity electrodynamics 
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The chaotic dynamics at the interaction between an ensemble of two-level atoms with an eigen- 
mode of a high-Q Fabry-Perot cavity and with an external amplitude-modulated field is considered. 
It is shown that in the case of an exact atom-field resonance and at the initial population of the 
atomic ground state, a Hamiltonian chaos in the system is transient. The difference in the chaotic 
dynamics for a Fabry-Perot and for a ring cavity geometries, as well as the conditions for the 
experimental observation of the transient Hamiltonian chaos are discussed. 



I. INTRODUCTION 



Tavis-Cummings (TC) model [1] which describes within the rotating-wave approximation (RWA) [2] the interaction 
between an ensemble of the two-level atoms with a high— Q cavity mode is a fundamental for quantum and nonlinear 
optics. Recently the basic dynamic effects predicted for the model have become available for the experimental 
verification (see [3,4] and references therein). 

Belobrov, Zaslavsky and Tartakovsky [5] showed that within the semiclassical approximation a breakdown of the 
^— I , RWA leads to the dynamical chaos (see also [6]). However, for the typical transitions in the optical and in the 
^ microwave domains, the violation of the RWA may occur only at rather high atom's density, when a gas model of the 
noninteracting two-level atoms is not more adequate. 

Nowadays other generalizations of the semiclassical TC-model with the Hamiltonian chaos in the RWA appeared 

^ [7-13]. These models describe the interaction of the three-level atoms with two cavity modes [7,8], and two-level 

atoms with one mode and with an external field injected into a cavity with a constant amplitude ^ [9-11], or with the 
: amplitude-modulated field [12]. For the last two cases, the interaction with an external field was taken into account 
■ within a standard spatially - homogeneous field approximation ("mean field model") [15]. 
^ , The common feature of all mentioned above models is the consideration of the sample shorter then the wavelength 
of the ring cavity with the only forward propagating wave. For the Fabry-Perot cavity, an account of spatial variation 
O , of the cavity field amplitude due to the interference between the forward and backforward propagating waves may 
Q ' influence sufficiently on the nonlinear characteristics of the atoms-field interaction. This fact has been well appreciated 
' in the theory of optical bistability (see [16, 17] and references therein). 

' In this paper we are concerned within a semiclassical approach the influence of the standing wave effects on the 
transition to the Hamiltonian chaos at the interaction between the two-level atoms with the field. We use the model 
[12], but instead of a ring cavity a Fabry-Perot cavity geometry is considered. 
It is shown, that in the simplest case of an exact atom-field resonance and at the initial population of the atomic 
^ ground state, the coupled Maxwell-Bloch system is reduced to the Hamiltonian system with 1.5 degrees of freedom — 
9^ . 8- periodically forced "Bessel pendulum" . Conditions for the transition to the dynamical chaos are found numerically. 
For the majority of the values of the external field amplitude and frequency, chaos is transient: the chaotic oscillations 
of the population difference and of the polarization become regular after some time. In contrast to other models of 
physical systems possessing the transient Hamiltonian chaos (see, e. g. [18 - 23]), in our case the effective potential 
of the Hamiltonian system with 1.5 degrees of freedom is of a multi-well type. This leads to rather complicated 
nonlinear dynamics of the system. In particular, a time interval of chaotic behavior has complicated dependence on 
the system's parameters. For trajectories with rather long chaotic part, a sensitivity of a character of settling of a 
regular oscillations of a cavity field with respect to small variations of the initial conditions or the parameters of the 
system is found. 

The paper is organized as follows. In Sec. 2 we introduce the model describing the dynamics of the system "two- 
level atoms -I- field" . The anzatz reducing the Maxwell-Bloch system to one equation of periodically driven "Bessel 
pendulum" in the case of the exact resonance and for some initial conditions is also considered in this Section. In 



^Note that this model may also describe the interaction of impurity center in a crystal with the light and phonons [f4]. In 
such case a cavity is absent. 
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Sec. 3 after a brief review on the trasicnt Hamiltonian chaos, we describe the nonhnear dynamics of our system. 
In conclusion, in Sec. 4 we compare the particularities in the dynamics of atoms-field interaction for a geometry 
of the Fabry-Perot and ring cavities, and discuss the possibilities of the experimental observation of the transient 
Hamiltonian chaos, and also point out some open problems. 

II. MODEL 

Our system consists of A'' identical two-level atoms with a transition frequency lvq and a dipole matrix element d, 

interacting with the radiation field mode with the frequency w « wq • The sample with atoms of the length L and the 
volume V completely fills a high— Q Fabry-Perot cavity {—L < z < 0). An external amplitude-modulated field E^xt 
is injected into the cavity at the plane z = 0. The field has the form 

Eext = EoF{t)cosLot, (1) 

where F{t) is a periodic function of time (| F{t) \= 1) slowly varying in comparison with the carrier frequency u>. 
In this paper we shall consider the simplest type of modulation F{t) = sin fit, but it can be expected that the main 
results are also applicable for the case, when the modulation F{t) consists of many harmonics. Following the classical 
paper of Spencer and Lamb [24], we shall consider the influence of an external field as a small local perturbation 
preserving a cavity mode structure of the form 

Esc{z,t) = Eit)smikz), (2) 

where k = n7rL~^ is the wave number of the chosen n-th cavity mode. The dynamics of atoms is described by the 
Bloch equations 

Si{z,t) = -L0t)S2{z,t), 

S2{z,t)=ojoSi{z,t) + ^Esc{z,t)S3{z,t), (3) 
2d 

S3{z,t) = -—Esc{z,t)S2{z,t), 

S! + Sl + Sl = l, (4) 

where the pseudospin variables Si {i = 1,2,3) are connected with the probability amplitudes for the population of 
upper ttj and lower bj levels of the j-th atom [9, 12, 25] 

Si{z,t) = ^^^i^v{0'*3^3 + ajb*), 
S2{z, t) = -±-Ef^^^{a*bj - ajb*), (5) 

S,{z,t) = ^^i:f^^y{\a^^\-\b]\). 

In (5) AV = (Az)7rr^ is a small averaging volume, z is the coordinate of slice center with the width Az ^ A (A is 
wavelength) , r is a characteristic radius of the sample with a gas of atoms, Ns is the number of atoms in the volume 

Ay (TV, » 1). 

The behavior of the field inside the cavity is governed by the Maxwell equation including an influence of the external 
fleld [24] 

E{t) + w^E{t) = -inP{t) - ^i^ext. (6) 



2 f° 

p{t) = - dzP{z,t)sm{kz), 

J-L 



where P{z, t) is a polarization of the two-level medium determined by 
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P{z,t) ^ ydSiiz,t). (7) 

The coupled Maxwell-Bloch system (3), (6), (7) completely describes the dynamics on the time interval shorter then 
all characteristic relaxation times of atoms and field. These equations may be simplified by separation of fast and 
slow dynamics. For this purpose make a transformation into the rotating system of coordinates 

Si{z, t) = u{z, t) cos(wt) — v{z, t) sin(wi), 

5*2(2, t) — u{z, t) sin(ajt) + t) cos{ujt), (8) 
Ssiz^t) = w{z,t), 

and introduce the envelopes of the self-consistent field in the following way 

E{t) = El (t) cos(wt) + E2 (t) sm{Lut) . (9) 
Assume that the conditions of a slow variance of the envelopes in a comparison with the carrier frequency u 

|u|<w|m|, \v\<^uj\v\, \w\<^oj\w\, 

I ^1,2 I ^1,2 I, (10) 

are valid. Then, substituting (1), (8) into (3), (6), (7) and using the RWA and a slowly varying envelope approximation 
(SVEA) [2] one can derive the following equations 

u — Av + w£2 sin(/cz), 
V = —Au + wSi sm(kz), 
w = ^{u£2 + vEi) si'n(kz) (11) 



2 r 

S^^ujI- dzv{z,t)sm{kz) + GF{t), 



oj2~ / dzu{z,t) sm{kz). 



In (11) A — uj — ujq is a detuning from the optical resonance, G = {ceo)/L, £j — dEj/h {j — 0,1,2), uJc — 
{2TTNd'^u>o/hVy/^ is a so-called cooperative frequency [26] characterizing the time of energy transfer between atoms 
and field in the absence of the external field (£^0 = 0)-O The Maxwell-Bloch equations possess the pseudospin length 
conservation law 

u{z,t)^ + v{z,tf +w{z,tf = 1. (12) 
Substituting (8) into (7), we obtain the representation for the polarization via the slow variables u and v 

Nd 

P{z, *) = — t) cos{ujt) - v{z, t) sm(u}t)\ . (13) 

The condition of applicability of the RWA and the SVEA may be written in the form 

max (^G^/^jWcf^^ < wq. (14) 

As a rule, for typical optical and microwave systems this condition is valid. 

Now the system (11) contains only slow variables. To solve it, it is necessary to define the initial spatial distributions 
of polarization components u{z,t = 0) = uo{z) and vo{z) as well as population difference ^0(2:) and field envelopes 
1^^1.2(0). Deriving a general solution of (11) looks as rather difficult problem, so we consider only particular solution 
for an exact resonance and at the specific initial condition. 



■^In modern literature also notions such as "collective Rabi frequency" and "vacuum- field Rabi frequency" are used for the 
determination of this frequency (see the discussion of the corresponding terminological questions in [4]). 
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It may be shown that for an exact resonance (A = 0), if u{z,t = 0) = £2{t = 0) = 0, then u{z,t) = £2{t) = for 
any t. ^ In this case it foUows from (12) that the polarization v{z,t) and the population difference w{z,t) may be 
parameterized by one angle variable (l>{z, t) 



Introduce the following anzatz 



v{z, t) = — sin0(z, i), 

w{z,t) = — cofi(l){z,t). (15) 

0(z,t) = i?(i)sin(fcz), (16) 
d{t) = f dt'£i{t'). 







From (15) one can derive 



v{z,t) = -2E^^oJ2„+i(i^(i))sin[(2n + l)fcz], 
w{z, t) = -Mm) - 2S^=o>^2„(t?(t)) cos(2nfcz), (17) 



where J„(a;) is the Bessel function. In these notastions the equations (11) may be transformed to one equation of the 
periodically driven "Bessel pendulum" 

^ + 2ujIJi{^) ^Gsinnt, (18) 
^ = £i =£{t). 

Taking into account that for small x: Ji{x) w x/2, one can see that the cooperative frequency lOc is just the frequency 
of small oscillations of the pendulum at G = 0. 

As far as we know, the anzatz (16) , (17) was first introduced in [27] when studying the metastable states in the 
coupled "atoms + field" system. The equation of the "Bessel pendulum" without the external driving (G = 0) has 
been also obtained in [27]. 

Now we discuss the problem of the limits of applicability of the anzatz (16), (17). Expansion (17) should be correct 
for any v(z, t) and w(z, t) including the initial distribution vo{z) and wq{z). It means that should exist such t?(0) = i?o 
, for which the expansions 

vo{z,t) = -2S,^^oJ2„+i(z9o)sin[(2n + l)fcz], 
wo{z,t) = -Jo(i?o) - 2I]^^oJ2„(z9o)cos(2nfcz), (19) 

are identities. ^ But the fulfillment of (19) is possible not for all fo(z) and wq{z). Among the initial distributions vq{z) 
and wo(z) for which (19) is correct, there are at least two physically significant cases: (i) the case of the periodic 
(5- like distribution and the population difference corresponding to z9o ^ 1, and (ii) the case of spatially homogeneous 
and weak initial excitation of two-level medium 'Wq{z) « — 1 and vq{z) k, 0. Such initial distribution corresponds to 
i?o ~ 0. Note that anzatz (16), (17) can be used for any initial values of the field f (0). 

In what follows, we shall consider physically most interesting case of atoms initially populated in the ground state 
{wq{z) = —1, wo(z) = 0), and when initially the field in the cavity is absent (f (0) — 0). It means that for the 
pendulum (18), the initial conditions are fixed: '!?(0) = ??(0) = 0. 

III. HAMILTONIAN DYNAMICS 

Equation (18) may be written in the Hamiltonan form 

dd _dH dp _ dH 
dr dp ' dr dd ^ 
H = //2 + V{d) - Gi? sin(OT), (20) 
V{d) = -2Jo{i9), 



^Or, if w(z,i = 0) = £i{t = 0) =0, then v{z,t) = £i{t) =0 for any t. 
*Prof. Wei-Mou Zheng called our attention to this fact. 
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where p = E/ojc, t = ujct, G = G/ojI, = il/ujc ■ The form of potential ^(1?) is shown in Fig. 1. 

Transfer into extended phase space [28] and introduce a pair of new canonicaUy conjugated variables {i/j, I) according 
to the formulas 

dAp_dH^ dl__aff 

dT~ dV dT~ d^}' ^ ' 

^ = TiT H = H + ni. 

The new Hamiltonian H is the integral of motion in the extended phase space {'d,p,il),I). We used formulas (21) to 
control the errors in the numerical calculations. 

In spite of rather simple form of the equation (18), the dynamics of this Hamiltonian system with 1.5 degrees of 
freedom may be very complicated. Such a behavior is due to the form of the potential ¥{■&) which is not periodic, and 
is a decreasing function of at | | — > (X) . Under the influence of a periodic forcing both regular and chaotic dynamics 
are possible. As V{'d) ^ at i? ^ icxo , then the asymptotic behavior of the dynamics is regular. Therefore, among 
with trapped periodic and chaotic trajectories, transient chaos is possible. 

In recent years the investigations on the transient Hamiltonian chaos attract considerable attention (see, c. g. , 
reviews [18,19] and references therein). The majority of papers deal with the study of particle scattering on 2D or 
2>D potentials [18-20]. It was quite recently realized that stochastic ionization of atoms and molecules or, generally 
speaking, any escape of chaotic trajectory from the potential well are also examples of the transient Hamiltonian 
chaos. The transient chaos under the conditions of stochastic ionization is not only much more poor studied, but 
also is organized more complicated then chaos in the potential scattering [22]. This is due to the fact that the phase 
space structure of the typical Hamiltonian systems is not uniform: Along with the stochastic seas the great number 
of islands with the stable dynamics are presented, which arc the source of the trajectory trapping for a long time. 

In paper [21] the stochastic ionization of kicked Morse oscillator is considered. Initial oscillator's state corresponds 
to the potential minimum. It was shown, that the time of chaotic escape from the well tesc may be sensitive to a 
small change of the parameters of the system. The boundaries of the regions with the different tesc are fractals in the 
space of parameters: amplitude — period of external driving. 

In paper [23] the escape of chaotic trajectory from 2D potential has been studied. It was shown that for strong chaos 
the directions and time of escape from the well are significantly sensitive to a small change of the initial conditions. 
The region in the space of the initial conditions corresponding to the trajectory escape during a fixed number of 
iterations appears to be a fractal. 

All mentioned above papers were devoted to the study of chaotic escape from 2D potential (2 degrees of freedom), 
or from ID potential which consists of one potential well (1.5 degrees of freedom). The particularity of our system 
with 1.5 degrees of freedom consists in multi-well structure of the effective potential (see Fig. 1). It should be noticed 
also, that in our case initial conditions are fixed just in the bottom of the first well. 

There may exist several types of the trajectories (both trapped and unbounded) : stable periodic, unstable periodic, 
trapped chaotic and unbounded chaotic demonstrating transient chaos. Due to nonlocal structure of our effective 
potential, the trajectory can be trapped not only within the first (central) well, but inside other wells too. Of course, 
time interval of trajectory trapping inside different wells is quite different. Trapped chaotic trajectories exist only for 
some singular values of parameters and the most typical behavior is transient chaos. During transient chaotic motion 
trajectory after escaping from the first well may randomly visit several wells and be trapped within these wells for 
some time. So, because of nonlocal character of the effective potential, transient chaotic motion may be organized 
rather complicated. 

We illustrate now the main types of the nonlinear dynamics. First, if the conditions of chaos are not satisfied, then 
a trajectory is located in the first well (—3.83 < d < 3.83) forever. A corresponding form of the field (momentum of 
effective system (20)) for regular motion is shown in Fig. 2. Under the conditions of chaos the behavior of the system 
is more various: 

(a) An example of unstable quasipcriodic trajectory trapped inside first (central) well during all time of computation 
is shown in Figs. 3, 4. Such trajectories are observed mainly for the parameter values belong the vicinity of the chaos 
boundary. 

(b) Examples of transient chaos are shown in Figs. 5 and 6. The trajectory quickly escapes from first well, visits 
randomly and "rotates" during some time inside few other wells and, by the end, a regular asymptotic motion sets 
in. (Note that minimums of the wells correspond to the values 1? w 0, ±7, ±13.3, ±19, ±25 etc.). The signs of the 
coordinate ?9 and the momentum p setting in at i — > 00 are random. Very small variation of the initial conditions 
leads to absolutely different asymptotic state (compare Figs. 5 and 6). This results is the interesting physical effect: 
the small fiuctuations of the initial cavity field produce the change of the asymptotic field state sign (Fig. 7). 
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(c) An example of chaotic trajectory trapped during all time of computation is shown in Fig. 8. The trajectory 
quickly escapes from the first well {tcsc ~ 10.8), then visits few neighbor wells and again returns to the first well. 
Corresponding chaotic field oscillations are shown in Fig. 9. 

(d) By the other hand, a transition to the regular behavior may be very quick and just after escaping from first well 
(Figs. 10, 11). 

The regions of regular and chaotic behavior in the space of parameters: amplitude of driving — frequency of 
modulation arc shown in Fig. 12. The chaos boundary is marked by stars, the regularity windows lying near this 
boundary are marked by boxes. The minimal value of the external field amplitude (G w 0.14) leading to chaos was 
observed at the frequency O ~ 0.85 . The trajectories with chaotic parts of different length exist in the regions 
situated above the chaos boundary. These trajectories have positive maximal Lyapunov exponent during the time 
interval of chaotic motion. Our preliminary numerical simulations show rather complicated dependence of the length 
of the chaotic phase on the parameters. But definitely the trajectories with long chaotic path (resc ~ 10^), as a rule, 
belong to the layer with the width ^0.1 near the boimdary with the chaotic component. 

We also studied the sensitivity of the sign of the self-consistent field asymptotic state with respect to the small 
variation of the external field amplitude. The results are presented in Fig. 13. Dimensionless perturbation parameter 
G was changing in the region [0.89 , 0.9] with a step 10~^. A regular asymptotic state with a positive field {£ > 0) is 
displayed by plus one ( +1 ), and the opposite case of the negative field {£ < 0) is marked by minus one ( -1 ). It is 
seen from Fig. 13, that in some regions of the perturbation parameter values G , the sign of the asymptotic field state 
is sensitive to small variation of G. Such kind of behavior is characteristic for trajectories with rather long chaotic 
path {Tesc Si 40-=- 100). 

In concluding part of this section we compare spatial structure of polarization and population difference at regular 

and chaotic dynamics (see formula (17)). An essential contribution into the expansion (17) make only spatial harmonics 
with the index of the Bessel functions less then its argument. For regular dynamics, the argument of Bessel function d 
is smooth regular function of time. In this case the spatial spectrum of v and w may contain many harmonics, but the 
energy transfer between different spatial modes is regular. By contrast, for chaotic temporal dynamics •& is random 
function of time. As a consequence, spatial chaos in the distribution of the population difference and polarization 
appears. 

IV. DISCUSSION 

It is interesting to compare a nonlinear dynamics of the system of two-level atoms + self-consistent field + external 
amplitude-modulated field for different geometry of the cavity. Generalized coupled Maxwell-Bloch system may be 
reduced to the periodically driven "Bessel pendulum" (18) in the case of the Fabry - Perot cavity. For the majority of 
paremeter values, Hamiltonian chaos is transient in this system. By contrast, for the ring cavity, the dynamics of the 
system within the same approximation is governed by the equation of the periodically driven physical pendulum [12] 

i^ + Lo'^siu'd = GF{t), 

where ^{t) is the Bloch angle. Due to periodicity of the corresponding potential V{"&) = — w^cost? and for periodic 
F{t), the dynamical chaos is stationary. Appearance of the transient chaos is possible only for a special choice of 
F{t), i. e. when F{t) contains zero harmonic [12]. Then, the effective potential loses its translation invariance. 

At present time the oscillatory energy transfer between atoms and field has been observed in both optical and 
microwave domains [4]. Present here some estimates. For dipole transition with d ~ IDebye and atom's density 
lO^^cm^"* , the value of the cooperative frequency ojc ^ lO^^s^^ may be large then the dissipative rates associated 
with the atom relaxation and cavity losses. We believe that the rapid progress in the experimental cavity quantum 
electrodynamics makes observation of the transient Hamiltonian chaos accessible. 

In summary, we considered the influence of the spatial structure of a standing wave on the chaotic dynamics in the 
interaction of an ensemble of two-level atoms with the field. It is shown, that for the majority of paremeter values, 
chaos is transient. For trajectories with long chaotic path, the asymptotic behavior of the field is sensitive to small 
variation of the initial conditions and parameters. Temporal chaos has manifestation in the spatial behavior of the 
polarization and of the population difference. 

In conclusion we point out some relevant problems for future research. First, a complicated nonlinear dynamics 
of periodically forced "Bessel pendulum" deserves more detail and deep study in the frameworks of general problem 
of transient chaos. Second, the problem of influence of quantum effects on chaos in potential scattering attracts 
great attention of researches (see review [29]). But the influence of quantum effects on the other types of the 
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transient Hamiltonian chaos is yet not considered. So, it could be interesting to study the influence of the cavity fleld 
quantalization on the transient chaos in the model discussed. 
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FIGURE CAPTIONS 

Figure 1. The form of the potential for the Bessel pendulum V{x) = —2Jq(x). 

Figure 2. Dependence of the self-consistent field p = £/wc on time for regular motion inside the first potential well: 
G = 0.8; n= 1.3. 

Figure 3. Weakly unstable chaotic orbit of length r = 200 located in the first well: G = 0.271; f2 = 0.8. 
Figure 4- Dependence of the field on time. The values of parameters are the same as in Fig. 3. 

Figure 5. Trajectory of the length t = 200 under the conditions of the transient chaos (G = 0.7; fl = 0.6). 

Figure 6. Same as Fig. 5, but with slightly changed initial conditions for the field: p(0) = 10""^. 

Figure 1. Changing of the asymptotic field state at small variation of the initial conditions. Process marked by 
boxes corresponds to the initial condition p(0) = 0, not marked process — p(0) = 10~^. Values of parameters are the 
same as in Figs. 5 and 6. 

Figure 8. Chaotic trajectory of length r = 200 at G = 0.6; ^ = 0.9. 

Figure 9. Chaotic field oscillations for the parameter values are the same as in Fig. 8. 

Figure 10. Quick escape from the first potential well (resc = 14.8) and setting of the regular motion at G = 0.4; 

n = 0.9. 

Figure 11. Dependence of the field on time. Values of parameters are the same as in Fig. 10. 

Figure 12. Regions of regular and chaotic motion in parameter space: amplitude — frequency. 

Figure 13. The influence of the variation of the driving amplitude on the sign of the asymptotic field state (for 
explanation, see main text), f2 = 0.4. 
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